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ABSTRACT 

We deal with the coupling of hyperbolic and parabolic sys- 
tems in a domain (1 divided into two disjoint subdomains Q + and 
Q ~ . Our main concern is to find out the proper interface conditions 
to be fulfilled at the surface separating the two domains. Next, we 
will use them in the numerical approximation of the problem. The 
justification of the interface conditions is based on a singular per- 
turbation analysis, that is, the hyperbolic system is rendered para- 
bolic by adding a small artificial "viscosity". As this goes to zero, 
the coupled parabolic-parabolic problem degenerates into the origi- 
nal one, yielding some conditions at the interface. These we take as 
interface conditions for the hyperbolic-parabolic problem. Actu- 
ally, we discuss two alternative sets of interface conditions accord- 
ing to whether the regularization procedure is variational or non- 
variational. We show how these conditions can be used in the 
frame of a numerical approximation to the given problem. Furth- 
ermore, we discuss a method of resolution which alternates the 
resolution of the hyperbolic problem within Ci~ and of the para- 
bolic one within fl + . The spectral collocation method is proposed, 
as an example of space discretization (different methods could be 
used as well); both explicit and implicit time-advancing schemes 
are considered. The present study is a preliminary step toward the 
analysis of the coupling between Euler and Navier-Stokes equa- 
tions for compressible flows. 1 


x This research was supported by the National Aeronautics and' Space Ad- 
ministration under NASA Contract No, NAS1-18107 while the second author 
was in residence at the Insitute for Computer Applications in Science and En- 
gineering (ICASE), NASA Langley Research Center, Hampton, VA 23665. 



1. INTRODUCTION 


In this work we deal with (initial) boundary value problems for partial 
differential equations (or systems) which change their character within the 
domain under consideration. Precisely, we consider problems which are of 
hyperbolic type in a subdomain Q~ of the whole domain Q and of parabolic 
type in the complement Q + . 

This interest is motivated by various applications. Among others, we 
emphasize the case of fluid dynamical problems for viscous, compressible flows 
in presence of a body, governed by the system of Navier-Stokes equations (see, 
e.g., [CM], [S]). A convenient numerical approach to the solution to these prob- 
lems relies upon the splitting of the physical region where flow occurs in two 
computational domains, one (say, Q + ) close to the body, where viscous terms 
are to be taken into account and another far away (say, Q~), where viscous 
terms may be neglected. 

This leads precisely to a coupled problem involving Euler equations 
(hyperbolic) in the far region and the complete Navier-Stokes equations in the 
near region. 

From a computational point of view, this splitting procedure carries obvi- 
ous advantages. In particular, we mention the possibility of using different 
solvers for the two subproblems. Of course, a crucial point in this framework is 
how to relate the two problems to each other at the interface separating the two 
subregions. This feature must be investigated for the differential problem, first: 
suitable conditions at the interface will be derived. Whatever numerical scheme 
is used, it must take these conditions into account. 

Identifying such conditions is generally understood whenever the two 
differential problems in the subregions are of the same kind. For instance, for 
the interaction of two second order elliptic problems, the interface conditions 
consist in requiring the continuity of the unknowns and of the flux (these con- 
ditions are transferred to any numerical scheme easily). Analogously, when 
coupling two Navier-Stokes problems, we must impose the continuity of the 
velocity and of the normal stress at the interface. Eventually, when coupling 
two hyperbolic systems of first order, we request the continuity of the unk- 
nown at the interface (unless the interface is a discontinuity line, in which case 
Rankine-Ugoniot equations ought to be fulfilled). 

When coupling Euler and Navier-Stokes equations, the proper interface 
conditions are not obvious, in advance. A possible way of deducing them is to 
see the coupled problem as a limit of two coupled Navier-Stokes problems with 
vanishing viscous terms in 0 - . 
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This approach can be adopted also in a simplified version of the problem, 
namely considering a coupling between hyperbolic and parabolic linear systems 
in one space variable (as well as their stationary counterpart). In this frame- 
work, the hyperbolic-parabolic problem is seen as limit of a coupling between 
two parabolic problems endowed with the usual transmission conditions at the 
interface. As we shall see, this procedure yields certain interface conditions for 
the limit problem. Although these conditions are somewhat physically meaning- 
ful, there is a sensible loss of continuity in passing from the approaching prob- 
lems to the limit one. In what follows we shall detail also a different type of 
limit procedure, which maintains a higher order of regularity for the limit prob- 
lem. 

To be more precise, let us state the problem we will discuss in this paper: 

In the interval (a ,c ), find w such that (h is any point internal to (a ,c )) 

w “ + A w~ + B w~ = F for x € = (a ,b ), t > 0, (1.1) 

w t + — (/t w x + ) x + A w + x + B w + = F for x € Q + = {b ,c ), t > 0, (1.2) 

with an initial condition and proper boundary conditions at x = a and x = c . 
Here A and B are two constant 3x3 matrices, while F is a given vector func- 
tion with three components (w is an unknown three dimensional vector); 
fi = fi{x ,t ) > /i 0 > 0 is a given viscosity. We assume that A has three real, 
non vanishing eigenvalues ( p 0 of them are positive and 3 —p 0 are negative): in 
particular, this implies that the system (1.1) is hyperbolic. 

For the above problem we are going to specify the interface conditions 
obtainable by the arguments previously mentioned. By the first approach, which 
we will refer to as variational , we find the following interface conditions at 
x = b , for all t > 0: 

T n w + = T n w - , (1.3) 

— /tiw x + + A w + = A w~: (1.4) 

the rectangular matrix T n has 3 — p 0 rows given by the left eigenvectors 
corresponding to the negative eigenvalues of A (see section 3.1). 

With the second approach, which we will refer to as nonvariational, the 
interface conditions at x = b and for all t >0 are: 

•w + = w _ , (1.5) 

Tn w x + ~ T n w x ". (1.6) 

In particular, note that (1.5) does imply continuity of all unknowns at the 
interface, while (1.4) gives continuity of the "flux" at the interface, allowing a 
discontinuity on the unknowns (actually, a mild discontinuity, as the jump has 
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the same order of the viscosity coefficient fx at the interface). 

The above results are presented in section 4, as a consequence of a pro- 
cedure of "increasing difficulty" carried out throughout sections 2 and 3. Pre- 
cisely, in section 3 we deal with the steady counterpart of (1.1), (1.2) and in 
section 2 we detail the coupling between two time independent equations, one of 
first order and the other of second order (the proofs of the abstract results are 
given in the Appendix). Although the problems of sections 2 and 3 might be 
regarded as autonomous problems, actually they are treated as intermediate 
steps toward the analysis of the main problem (1.1), (1.2). For each and every 
problem, we present the numerical approximation based on the spectral colloca- 
tion method and show how the interface conditions are used in this frame. This 
could be done for numerical methods based on different approaches, as well. 
Here we just remark that, in the numerical scheme, we must supplement the 
above interface conditions suitable compatibility relations at the interface. These 
arise from the hyperbolic nature of the problem in Q~: a thorough discussion is 
made in sections 2.2, 3.2, 4.2. 

We end this introduction by noticing that (1.1), (1.2) present some simi- 
larities with the coupling between Euler and Navier-Stokes equations we men- 
tioned at the beginning as a driving motivation for our work. The relevant 
difference lies in that the viscous terms in Navier-Stokes equations do not enjoy 
the particular diagonal structure as in the right hand side of (1.2). Since our 
analysis relies heavily upon this feature, there is no immediate application of 
our results to the coupling between Euler and Navier-Stokes equations. 
Nevertheless, it seems that several elements of our approach can be useful in 
that problem, too. From this point of view, the present work is an intermediate 
step toward our goal. 


2. HYPERBOLIC-ELLIPTIC INTERACTION: THE SCALAR CASE 

In this section we consider a one dimensional, linear, scalar problem. The 
two subsections are devoted to the analysis of the continuous problem (with 
special concern to different elliptic regularizations) and to its numerical approxi- 
mation, respectively. 


2.1. The differential problem 


We begin by stating the boundary value problem, as follows. Let 
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(i) a, b , c be real numbers, with a < b < c ; 

(ii) <x,)3,/z be functions defined in [a,c], with a ^ 0, /u.(x )^yu 0 > 0 for 
x € [a ,c ]; 

(iii) / be a function defined in [a ,c ]. 

Then, consider the problem 

(P): find u defined in [a ,b ], v defined in [6 ,c ] such that 


a- u x + /3 u = / 

in (a ,b); 

(2.1) 

(/zv x ) x + a v x + p v = / 

in (b,c); 

(2.2) 

v (c ) = 0; 


(2.3) 

u (a ) = 0, 

if a > 0 in [a ,c ]. 

(2.4) 


Clearly, the formulation of problem (P) is incomplete: it needs one coupling 
condition between u and v at the interface b , when a > 0 in [a ,c], while two 
coupling conditions are required if oc < 0 in [a ,c] (in this case, (2.4) does not 
hold). Moreover, we may allow (2.3) to be substituted by v x (c ) = 0, if a > 0 
in [a ,c ]. 

Remark 2.1 Problem (P) may be regarded as a stationary problem (in this case 
3 might vanish identically) or else as a time discretization of an evolution 
advection-diffusion problem (hyperbolic in (a ,b ) and parabolic in ( b ,c )) by 
an implicit method (in this case, /3 behaves essentially like the reciprocal of 
the time discretization step). For this reason, we will always refer to prob- 
lem (P) as to a "hyperbolic-elliptic" problem, even if (P) is a purely steady 
problem. By the way, we just note that the characteristic lines of the evo- 
lution hyperbolic problem enter the domain (a ,b ) X (0, + oo ) across 
{a } x (0, + oo ), when a > 0 and across {b } x (0, + oo ), when a < 0. This 
is the reason why we choose to impose condition (2.4) among others, which 
are equally admissible for the time-independent problem. When or < 0, the 
same argument suggests not to impose any boundary condition at x = a 
(though admissible for the very equation (2.1 )); on the contrary, we are led 
to consider a condition on u at x = b. In the frame of the global problem 
(2.1), (2.2), this condition rea’ds as an interface condition. 

Two different types of elliptic regularizations are possible for problem (P), both 
acceptable for some reason. We will see that the two ways are essentially 
different as for the behavior at the interface. 
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The case a > 0. 


Given e > 0, consider the problem 

(P € ): to find u e defined in [a ,b ], 
- eu exx + a u €iX + 3 u e = f 
- (MV e , x ) x + oc v ex + 3v e = / 

u € (a ) = 0; 
v € (c ) = 0; 

( i ) u € = v € 

(“ ) eu e,x = P- v e,x 


: defined in [b ,c ] such that 


in (a ,b ); 

(2.5) 

in ( b,c ); 

(2.6) 


(2.7) 


(2.8) 

at x = b. 

(2.9) 


(P 6 ) is equivalent to a variational problem on the whole of (a ,c ); condition 
(2.9) expresses that u e and v e join continuously at b and that the flux across b 
is continuous, too. 

About the existence of solutions to problem (P e ) and their behavior as 
t — ► 0, the following result holds (see Appendix, where the appropriate choices 
of functional spaces are made and the regularity assumptions on the data are 
specified). 


Proposition 2.1 Assume the coerciveness condition 

2 3 — oi x > 0 in [a ,c ]. 

Then, problem (P e ) has a unique solution. Furthermore , as 
converge to a pair of functions u , v which satisfy (2.1), 
and the interface condition 

a ' u =■ — p.v x + oi v at x = b. 


( 2 . 10 ) 

0, u e and v e 
(2.2), (2.3), (2.4) 

( 2 . 11 ) 

□ 


Remark 2.2 (2.1 1) means that the flux across b is conserved, as e-* 0. On the 
contrary, analytical solution of (P 6 ) shows that u and v do not join con- 
tinuously at b , in general. Actually, the closed form of the solution (as 
well as numerical experiments, see subsection 2.2.3) shows that the jump 
between u and v at b has the same order as /i, when p. — ► 0. 
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A second approach is to consider the following problem 

(Q € ): to find u e defined in [a ,£>], v € defined in [ b ,c] such that (2.5), 
(2.6), (2.7), (2.8) and (2.9i) hold, along with the condition 

u €x = v ex at x = b. (2.12) 

(Oe) is equivalent to a nonvariational elliptic problem on the whole of (a ,c ): 
now we look for a pair of functions u e , v £ which have a C 1 junction at b . 

Proposition 2.2 Assume the coerciveness condition (2.10). Let u € , v e solve 
problem (Q e ). As e — > 0, u e and v € converge to a pair of functions u , v 
■which satisfy (2.1), (2.2), (2.3), (2.4) and the continuity condition 

u(b)=v(b) (2.13) 

at the interface. 

□ 


We remark that (2.12) is not preserved, in general, as e— ► 0: this can be checked 
on the closed form of the solutions to problem (Q e ), in some particular cases. 
Moreover, this feature is clearly shown by the numerical results presented in 
subsection 2.2.3. Thus, we are approaching a solution to problem (P) which is 
continuous but not C 1 at b . 


The case a < 0. 


In this case, one can consider the same problems (P € ) and (Q e ) as before. 
However, for a reason which will be clear in section 3, we prefer to perform a 
slight change in the two problems, namely replacing the Dirichlet condition 
(2.7) with a Neumann one. Note that the original problem (P) has no condition 
at all for x = a . Thus, we are dealing with a new couple of problems, which we 
denote by (P e ) N and (Q € ) v . respectively. For clarity, we state them in detail. 

(P £ )a^ : to find u € definal in [a ,b ], v € defined in [b ,c ] such that (2.5), 
(2.6), (2.8) and (2.9) hold, along with the condition 

u e,x( a ) = 0 . 


(2.14) 
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(OJ.V : to find u € defined in [a ,b ], v £ defined in [b ,c] such that (2.5), 
(2.6), (2.8), (2.9i), (2.12) and (2.14) hold. 

The difference with respect to the case a > 0 lies in the asymptotic behavior 
and, more precisely, in the interface conditions (remind that the limit problem 
(P) needs two conditions at b , in this case). The abstract analysis shown in the 
Appendix yields the following results (again, we do not specify the regularity 
on the data and on the unknowns here). 

Proposition 2.3 Assume the coerciveness condition (2.10). Then, problem (P € )/y 
has a unique solution. Furthermore, as e — ► 0, u € and v € converge to a pair of 
functions u , v which satisfy (2.1), (2.2), (2.3) and the following interface 
conditions: 

u(b)=v(b), (2.15) 

v x (b) = 0. (2.16) 

□ 


Proposition 2.4 Assume the coerciveness condition (2.10); moreover, suppose 
that j3 > 3o> 0 in l<3 ,^> ]. Let u € , v € solve problem (0 6 )^- As €-> 0, u e 
and v € converge to a pair of functions u , v which satisfy (2.1), (2.2), (2.3) 
and the following interface conditions: 

u(b)=v(b), (2.17) 

u x (b) — v x (b ). (2.18) 

□ 

We point out that the condition at a for both (P € ) A > and (Q e )^ is lost in the 
limit, as it is natural for this kind of problems. 

Remark 2.3 By means of both approaches, the limit functions u and v enjoy a 
continuous junction at b . But the derivatives behave in a very different way 
(see (2.16) and (2.18)). Indeed, the limit of the solution to (P € )a? shows an 
angle at b , in general, while the limit of the solution to (Q € ) A is C 1 at b. 
Thus, as in the previous case, the nonvariational approach is able to preserve 
an order of regularity higher by one, with respect to the variational one. 


Remark 2.4 The two regularized problems with the original Dirichlet condition 
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(2.7) have the same type of asymptotic behavior as the problems with the 
Neumann condition (2.14). The difference lies in that in the Dirichlet case 
the value u e (a ) does not converge to the corresponding value u (a ), which is 
true for the Neumann case of problems (P e )jv and (0«^ • 

Remark 2.5 A comment is needed about (2.18). This condition calls into play 
the first derivative of the solution to (2.1) at b: but (2.1) is a first order 
equation, hence (2.18) involves a boundary operator of the same order as the 
interior equation. Thus, the left hand side of (2.18) must be compatible 
with the collocation of the equation (2.1) at b . Precisely, whenever the data 
are smooth, we expect equation (2.1) to hold at b, hence (2.17) and (2.18) 
imply 

civ, + |3 v = / at x = b. (2.19) 


2.2. The numerical approximation 

Set Q - = ]a,M, Q + = On the reference interval [—1,1], let us con- 

sider the Chebyshev collocation points 

x*=- cos , j = 0, • • • ,N , (2.20) 

whose images in the interval 11 x are denoted by [x * }. Note that Xq = a, 
X N = x 0 = b , Xjf = c . 

As an initial step, we consider two separate boundary value problems: a 
first order problem in fl” and a second order elliptic problem in ft + . Next, we 
introduce their numerical approximations based on the spectral collocation 
method. This presentation has the aim of providing the reader a guideline to the 
numerical approach of the coupled problem (P). 

2.2.1. The split model problem 

The two separate differential problems in ft - and ft + are the following 
(we keep the same terminology as in section 2.1). 

"Hyperbolic” problem in Q~: 

au x + fru = f i n Q- t 

u(a)= u a if a > 0, (2.21) 

u{b) = u b if oc < 0, 
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where u a and u b are given. The motivation for the different choice of boundary 
conditions is given in Remark 2.1. 

Elliptic problem in Q + : 

- ) x + a v* + 0 v = f i n n+, 

B b v = v b at x = b , (2.22) 

B c v = v c at x = c , 

where v b and v c are given and B b v and B c v are suitable combinations of v and 
v x leading to a well posed problem. 

The spectral collocation approximation to (2.21) is as follows (see, e.g., 
[CHQZ], Ch. 10 and 11). We look for u N € (the space of algebraic polyno- 
mials of degree ^ N ) such that 


a u N ,j+0% = / 


at Xj , j = 1, • * • ,N — 1, (2.23) 


supplemented by the two boundary equations: 


if oc > 0 


if a < 0 


(i ) u N = u a at x 0 , 

Ui) (X Ufy x f} lift — f Ctt Xfl t 

(i ) oi u N x u N = f at x 0 , 
(ii ) ujy = u b at xff. 


(2.24) 


(2.25) 


The numerical approximation to (2.22), based on the spectral collocation 
method, is as follows. We look for v N € P N satisfying 


[ I n ip- v N ,x + <* V N ,x + $ V N — f 

at Xj + , j = 1, • • • ,N — 1, 

(2.26) 

B b V N = v b 

at Xq , 

(2.27) 

& 

•e 

II 

at xj , 

(2.28) 


where 1 N is the interpolation operator at the points x+. 
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2.2.2. The original coupled problem 


Now we are in a position to describe the numerical approximation to the 
original coupled problem (P), taking (2.23 )-(2 .28 ) into account. 

1. At the interior points of Q~ and Q + , we impose the set of equations 
(2.23) and (2.26), respectively. 

2. At x = a , we impose either (2.24i) (with u a = 0) or (2.25i), according to 
the sign of or. 

3. At x = c , we always enforce v N = 0 (which corresponds to (2.28) with 
v c = 0 and B c v N = v v ). 

4. At x = b , we need two equations, in order to close the algebraic system. 
These depend both on the sign of a and on the interface conditions pro- 
vided by either elliptic regularization (see section 2.1 ). In particular: 

(a) if a > 0, we impose (2.24ii), along with either 

— pv N x + a- v N = o u N ( variational approach ) (2.29) 

or 

v n ~ u N ( nonvariational approach ); (2.30) 

(b) if a < 0, we impose the condition 

= (2.31) 

(i.e. (2.25ii), with Uf, = v A -(x 0 + )); the remaining equation is given by 
either 

Vftx = 0 ( variational approach ) (2.32) 

or 

v A’,x = U N ,x ( nonvariational approach ). (2.33) 


We note that (2.29), (2.30), (2.32) and (2.33) are but special versions of (2.27), 
with suitable 'choices of B b and v b . These are specified in table 1, which sum- 
marizes the equations to be fulfilled by the numerical solution at each colloca- 
tion point (including boundary and interface). 
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Collocation 

points 

j cx>0 

1 

1 

O 

V 

a 

variational 

nonvariational 

variational 

1 

nonvariational j 

x = a 

(2.241) 

(2.24i) 

(2.23) 



1 

(2.23) 

xpK 1 </V-l 

(2.23) 

(2.23) 

(2.23) 

(2.23) 

x = b 

(2.27) (1) (L) 
(2.23) (C) 

(2.2 7 ) (2) (L) 
(2.23) (C) 

(2.25U) (L) 
( 2.2 7 ) (3) (L) 

(2.25ii) (L); 
(2.27) U) (L) 

i 

V/ 

V/ 

(2.26) 

(2.26) 

; 1 

(2.26) 

i 

(2.26) 

r 

•X = c 

(2.28) 

i 

(2.28) (2.28) 

(2.28) 


Table 1. Numerical approximation to problem P by spectral collocation method : (L) = 
limit condition given by the asymptotic analysis ; (C) = compatibility condition 
(the transport equation must be collocated at the outflow boundary of flj ); 

— fx Vtf x + av N and v b — ocu N (x^ (1) ) ; 

and v 6 = u N (x r j l) ) : 

Vtf tX and v b = 0; 
v N.x and v b = u N x (x f } 1) ) . 


(1) 

with 

B-v N 

(2) 

with 

B-V N 

(3) 

with 

B~v n 

(4) 

with 

B-v n 
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2.2.3. Some numerical results 


Now we present several numerical experiments which support the theoreti- 
cal results obtained in the previous subsections. We deal with the elliptic regu- 
larizations of problem (P), taking (a,c)= (—1,1), with b = 0. In all cases 
(<* > 0 or a < 0, variational or non variational approach), the equations are 


- + a Ue,x + 3 U e = f 


in (—1,0); 

~ (/ iV €,x )x + a v €x +|3v £ = / in (0,1). 

The interface conditions change according to the regularization chosen: 


{variational ) 


(i ) = 

( H ) = fi v 


at x = 0 


€,X 


(2.34) 

(2.35) 


(2.36) 


or 


( nonvariational ) 


U) 


u e = v. 


(« ) u e,x = 


at x — 0. 


(2.37) 


€ ,X 


The boundary conditions will be distinguished later. 


These problems are solved by the Chebyshev collocation method described 
in advance for fully elliptic problems of the form (2.22). 

(To be more precise, we have implemented the collocation method in a 
domain decomposition framework, in order to achieve the highest precision. To 
this end, three subdomains are used; within each of them, we take 50 points; the 
middle subdomain includes the interface point x = 0. At each interface between 
subdomains the C 1 continuity is enforced directly (see [FQZ]).) 


The data we have used are the f ollowing 


a = —1, b = 0, c = 1, / = 1, ii = 1, j3 = 1 +x 2 . (2.38) 


A homogeneous Dirichlet condition is enforced at x = 1. About the point 
x = — 1, we consider the case of a homogeneous Dirichlet condition, to begin 
with. 


In Figure 2.1 we graph the results obtained for the variational approach, 
with e= 0.005 and e— 0.1, when a = 1. In agreement with our theoretical 
results (see Remark 2.2), the solution exhibits a discontinuity as e— » 0 at the 
interface point x = 0. The discontinuity is revealed by the presence of oscilla- 
tions near the interface, due to the Gibbs phenomenon. However, the jump is of 
the same order as the viscosity coefficient fi, as shown in Figure 2.2. 

Figure 2.3 displays the results obtained for the nonvariational case, using 
the same data as in Figure 2.1. Note that, as 0, the solution is continuous 
(though not C 1 ) at the interface point, as predicted by (2.1 7). 
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The comparison between variational and nonvariational approaches is 
clearer in Figure 2.4, where we take e = 0.005. 

In Figures 2.5 and 2.6 we present the results obtained using the two 
approaches, with the same data as before, but with a = — 1. As predicted by the 
theory (see Propositions 2.2 and 2.3), as € — > 0 the nonvariational solution 
remains C 1 , while the variational one is just C°. 

Finally, Figure 2.7 reports the results obtained with <x = — 1, € = 0.005 
and a homogeneous Neumann condition at x = — 1 (rather than the Dirichlet 
one), using the variational and nonvariational approaches. 


Figure captions 


Figure 2.1 Results for the variational approach , with a = 1: €=0.1 (dashed 
line), € = 0.005 (solid line). 

Figure 2.2 Results for the variational approach , with or = 1 , €= 0.005/ 
p — 1 (solid line), p = 0.1 (dash-dot line), p = 0.01 (dash-dash 

line) 


Figure 2.3 Results for the nonvariational approach, with a = 1: €=0.1 

(dashed line), e= 0.005 (solid line). 

Figure 2.4 Comparison between the two approaches, with a= 1, €= 0.005: 
nonvariational approach (dashed line), variational approach (solid 

line). 

Figure 2.5 Results for the variational approach, with oe = — 1: €=0.1 

(dashed line), e= 0.005 (solid line). 


Figure 2.6 Results for the nonvariational approach, with a = — 1: €=0.1 
(dashed line), e= 0.005 (solid line). 


Figure 2.7 Comparison between the two approaches, with a = — 1, €= 0.005, 
with homogeneous Neumann boundary condition at x = — 1 arid 
homogeneous Dirichlet boundary condition at x = 1: nonvaria- 
tional approach (dashed line), variational approach (solid line). 






Fig. 2.3 
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Fig. 2.4 
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Fig. 2.5 
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3. HYPERBOLIC-ELLIPTIC INTERACTION: THE (TIME INDEPENDENT) 
VECTOR CASE 

In this section we consider a boundary value problem for a system of three 
linear equations. Precisely, we deal with the stationary problem associated to 
( 1 . 1 ), ( 1 . 2 ). 


3.1. The differential problem 

With the notations of the introduction, we seek for a pair of three dimen- 
sional vector f unctions w - and w + such that 

A w “ + B w~ = F in fl - , (3.1) 

— (/zw/)* + A w x + + B w + = F in Q + . (3.2) 

About boundary conditions, we must distinguish between the points a and c . 

At x = a , we prescribe exactly p 0 conditions on w~, where po is the 
number of positive eigenvalues of the matrix A . These conditions take the form 

G“w“ = q - at x = a , (3.3) 

where q _ is a given vector with p 0 components, while G~ is a po* 3 matrix 
with rank p 0 . The choice of G~ is subject to some restrictions that will be 
specified later. 

At x = c , the boundary conditions can be written in the general form 

G + w + + H + w x + = q + at x = c , (3.4) 

where G + and H + are 3 x 3 matrices and q + is a given vector with three com- 
ponents. (3.4) must provide 3 independent equations, which are admissible for 
the elliptic system (3.2). In general, (3.4) yields a coupling between the three 
components of w + and their derivatives. However, in some special cir- 
cumstances, (3.4) might lead to three equations, each of them containing only 
one component and/or its derivative. 

Problem (3.1 )-(3.4) needs (6 — p 0 ) further conditions at the interface point 
x = b. Essentially, three of them are requested by the elliptic system (3.2), 
while (3 — po) (the number of negative eigenvalues of A ) pertain to the hyper- 
bolic system (3.1). 

To write down the interface conditions, let us introduce the matrix T 
which diagonalizes A and denote by 


f 
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A = T A r~ l 


(3.5) 


the (diagonal) eigenvalue matrix. We write A as 


A = 


A p 0 
0 A„ 


(3.6) 


where A p is the diagonal matrix of the p 0 positive eigenvalues of A , while A n 
is made by the remaining 3 — p 0 negative eigenvalues. Correspondingly, we 
write T as 


T = 


T 

L p 

T 

1 n 


(3.7) 


where T p is the submatrix of the first p 0 rows of T and T n is the rest (note 
that the rows of T are made by the left eigenvectors of A ). 


The interface conditions we consider here are of two types: either 


(i ) — /iw/ + A w + = A w (3 cond.s ) 

(ii ) T n w + = T n w _ (3 — p o cond.s ) 


( variational approach) or 

(i ) w + = w~ (3 cond.s ) 

(ii ) T n w/ = T n w x ~ (3 — p 0 cond.s) 

(nonvariational approach). 

In both cases, we impose as many conditions as requested. 

The rest of this subsection is devoted to a mathematical justification of 
(3.8) and (3.9) by means of the asymptotic procedure on elliptic regularizations, 
in analogy with the scalar case. 

Precisely, for a given e > 0, we consider the regularized problem 


— €w“ xx + A w“ x + B w“ = F in fi _ , (3.10) 

— (p w € ^ x ) x + Aw € + x + B w € + = F in ft + , (3.11) 

with boundary conditions 

G"w“ = q" 

„ at x = a , (3.12) 


G + w € + + H + w e ^ x = q + at x = c 
and interface conditions at the point x = b 


(3. i 3) 
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or 


( variational 
approach ) 


w £ = w e + 

€ W “ x = /X w+ x 


(3.14) 


{non variational 
approach ) 


w e = w € + 
W C Tx = 


(3.15) 


In (3.12), the original boundary condition has been added a homogeneous Neu- 
mann condition on T n w e ~: this is not the only possibility, but it is optimal, in 
some sense (see Remark 3.1). 

In order to exploit the results of section 2.1, it is natural to diagonalize the 
system (3.10), (3.11). This is done by introducing the characteristic variables 
associated to the system, namely 


Zz = T w € * 


denote by (z/ ) p the first p 0 components of z* and by (z € * )„ the remaining 
3 — p 0 components. Thus, (3.10M3.15) imply that z e * satisfy the equations: 


“ ez € xx + Az €> x + B t z € = F r in Q , 

“ (M z« + x )x + Az+ X + B t z € + = F r in Q + 

(where B T — T B T -1 and F T — T F), boundary conditions 


(3.16) 

(3.17) 


G r z e = q 

( - s „ at x = a , (3.18) 

(Z € ,x )n = 0 

Gf Ze + + Hjz£ x = q + at x = c (3.19) 

(where G/-* = G* T -1 and /7 r + = H + T ~ 1 ) and interface conditions at the 
point x =■ b 


{variational 
approach ) 


z € = z + 

fZ-lx - fJ-Z-Zx 


(3.20) 


or 


{nonvariational 
approach ) 


z e = z e + 


-*€,X 


= Z, 


€,* * 


(3.21) 
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In (3.18), the matrix Gy must satisfy the following assumption: 

the submatrix given by the first (3 22) 

Po columns of Gf is nonsingular. 

(3.22) poses restrictions on the choice of G - in (3.3), depending on A . 

We are now in a position to use Propositions 2. 1-2.4, whence we get the 
following convergence results (we suppose existence of solutions to the regular- 
ized problems). 


Proposition 3.1 (Variational approach) As 0, the solution z e , 

(3.20) converges to a pair of functions z~ , z + which satisfy: 

, z € + to (3.16)- 

Az x + B t z — 

Fy in , 

(3.23) 

~ (A z x + ) x + Az x + B t z + = 

Fy in Q, + , 

(3.24) 

O 

N , 

II 

1 

at x = a, 

(3.25) 

Gf z + + H T + z + = q + 

& 

x 

II 

o 

(3.26) 

M z p x + A p z p A p z p 

at x = b , 

(3.27) 

Zn + = Z " 

II 

Cr 

(3.28) 

O 

II 

« 

.O 

II 

X 

13 

(3.29) 


(in (3.27)-(3.29), z p ± denotes the first p 0 components of z ± and zf denotes 
the remaining 3 — p 0 components). 

□ 


Proposition 3.2 (Nonvariational approach) As e— » 0, the solution z e , z e + to 
(3.I6)-(3.19) and (3.21) converges to a pair of functions z~, z + which satisfy 
(3.23)43.26) and 

z + = z~ at x = b , (3.30) 

Z n + ,* = Zn~* at X = b. (3.31) 

□ 

By re-transforming these results in terms of the physical variables w* , we find 
(3.1)-(3.4) and either 

(i) - fxT p v/? + A p T p vs + = A p T p w~ ( p 0 cond.s ) 

(ii ) T n w + = T n w~ (3 — po cond.s ) (3.32) 

(Hi ) T n w x + =0 (3 — p o cond.s ) 
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or (3.9), according to the regularization chosen. It remains to show that (3.32) 
is equivalent to (3.8). Actually, taking (3.32ii) into account, (3.32iii) can be 
written as 

- V-T n w x + + A n T n w + = A„ T n w". 

Together with (3.32i), this last condition gives 

— fiT w x + + AT w + = AT w~. 

Multiplying by T~ } and recalling that p. is a scalar function, we get 

— /zw x + + A w + — A w _ , 

whence (3.8) follows. 

Remark 3.1 A Dirichlet condition on T n w~ in (3.12) (i.e. a Dirichlet condition 
on (z € ") n in (3.18)) is as good as the Neumann condition we considered in 
(3.12), provided (3.25) involves the p 0 characteristic variables correspond- 
ing to positive eigenvalues only. This means that the last 3 — p 0 columns of 
Gf ought to vanish identically. Essentially, the reason of this drawback is 
that the Dirichlet condition cannot guarantee the convergence of (z<T) n (a ) 
to z ). Thus, the strategy of reducing the analysis of the system to that 
of the scalar case cannot deal with a condition of type (3.3) involving the 
value of ( a ). However, a more sophisticated vector approach could be 

performed, capable of overcoming this difficulty (see [BBB], [BR], [L]). 


3.2. The numerical approximation 

We adopt the notations of section 2.2 for the collocation points. 

The spectral collocation approximation to problem (3.1)-(3.4) reads as fol- 
lows. We look for € (P^ ) 3 and € (P/y) 3 satisfying: 

A w.y x + B = F at x~ , j = 1, • • • ,N— 1, (3.33) 

— [/ A r(/iw^ x )] x -l- A w^ >x + B = F at x + , j = 1, • • • ,N — 1. (3.34) 

The conditions at x = a are of two types: 

(i) p o prescribed boundary conditions (see (3.3)): 


G w iV = q 


at x 0 , 


(3.35) 
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(ii) (3 — Pq) compatibility conditions: 

T n [A ' X + B — F ] = 0 at Xq . (3.36) 

Note that (3.36) are nothing but the collocation at Xq of the equations on the 
characteristic variables corresponding to negative eigenvalues: they generalize to 
the vector case the compatibility condition (2.25i) for the scalar case, yielding a 
stable and consistent scheme (see, e.g., [CQ]). 

At the right boundary point c, we enforce the prescribed boundary condi- 
tions (3.4) on the discrete solution, namely 

G + Wyy + x = q + at x£ . (3.37) 

Now, we come to the conditions at the interface point b . As usual, we dis- 
tinguish between the variational and the non variational approaches which have 
been used. The results of the analysis presented in section 3.1 (see Propositions 
3.1 and 3.2) suggest the proper continuity conditions to be enforced at the inter- 
face point. 


(a) Variational approach. 

(i) p o compatibility conditions on the equations corresponding to the positive 
characteristic variables: 

T p [A x + B w^ — F] = 0 at b (= x £ ); (3.38) 

(ii) (3 — Pq) conditions of continuity on the negative characteristic variables: 

T n w,v = T n Wjv at b (= x ), (3.39) 

obtainable from (3.8ii); 

(iii) 3 conditions of continuity of the "flux" on the physical variables: 

+ Awjy = at b i= Xq ). (3.40) 

obtainable from (3.8i). 

Remark 3.2 Notice that the hyperbolic system (3.33) has been supplemented 
three conditions at the interface point b (= x N ) in (i) and (ii) (see (3.38) 
. and (3.39)). Similarly, the elliptic system (3.34) has been given three 
Newton-like conditions at the interface point b (= Xq ) in (iii) (see (3.40)). 
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(b) Nonvariational approach. 

Ci) p o compatibility conditions on the equations corresponding to the positive 
characteristic variables, given by (3.38); 

(ii) (3 — p 0 ) conditions of continuity on the negative characteristic variables, 
given by (3.39); 

(iii) p 0 conditions of continuity on the positive characteristic variables: 

~ r p ™N at b (= *o + ): (3.41) 

both (ii) and (iii) are obtainable from (3.9i); 

(iv) (3 — p 0 ) conditions of continuity of first derivatives on the negative 
characteristic variables: 

T n w — T n Wy x ax b (= Xq ), (3.42) 

obtainable from (3.9ii). 

The same kind of considerations as in Remark 3.2 can be made in this case, too. 
We note that (ii) and (iii) amount to require that Wjv = at b. 

Remark 3.3 An efficient (and quite natural) method to solve problems of the 
form (3.33M3.37), supplemented with the interface conditions (3.38)- 
(3.40) (or (3.38), (3.39), (3.41) and (3.42)), relies upon an iterative pro- 
cedure alternating the solution of a hyperbolic problem in Q~ and of an 
elliptic one in ft + . 

At each step, the iterative method entails within Q~ the solution of the 
hyperbolic problem (3.33) with the boundary conditions (3.35) and (3.36) 
at the left hand boundary Xq , and (3.38), (3.39) at the right hand boun- 
dary Xfl . Next, in Q + we solve the elliptic problem (3.34) with the boun- 
dary condition (3.37) at the right hand boundary x and the conditions 
(3.40) (or (3.41), (3.42)) at the left hand boundary Xq . Finally, a relaxa- 
tion procedure on the interface variables is generally needed, in order to 
ensure the convergence of the above process. 

The details and the convergence analysis will be presented in a forth- 
coming paper. 
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4. HYPERBOLIC-PARABOLIC SYSTEMS FOR TIME DEPENDENT PROB- 
LEMS 

In this section we consider the problem (1.1), (1.2) presented in the intro- 
duction, endowed with its boundary, initial and interface conditions. 


4.1. The differential problem 

With a , b , c chosen in the usual way, we look for a three dimensional vec- 
tor valued function w* defined for x € Q* , r >0, satisfying 

Wj“ + A w “ + B w _ = F for x € O - = (a ,b ), t > 0, (4.1) 

w, + — (yuw x + )* + A w x + + B w + = F for* € 0 + = (b ,c ), t > 0, (4.2) 

where A , B , F and yu are given as like as in the introduction. 

The system (4.1), (4.2) must be given an initial condition 

w ± (;t ,0) = (;r ), *€0* (4.3) 

and boundary conditions, which we take again of the form (3.3) and (3.4), 
namely 

G~ w - = q“ at x = a , t > 0, (4.4) 

G + w + + H + w x + = q + at x — c ,t > 0, (4.5) 

where G~ , G + , H + , q~ and q + may depend on t . 

Analogously, at the interface line {b}x(0,+oo) we impose conditions 
which are the natural extension of (3.8) and (3.9) to the evolution case: either 


( variational 
approach ) 


Cl) T n w + = T n w 

(ii ) — /iw, + + A w + = A w - 


(4.6) 


or 


( nonvariational 
approach ) 


(z ) w + = w 
(«■) ^n w x + = r n w; 


(4.7) 


for x = b and for t > 0. 

The interface conditions (4.6) or (4.7) might be derived directly by means 
of regularized parabolic problems, in analogy to the procedure presented in sec- 
tion 3.1. 
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On the other hand, several heuristic justifications of these conditions may 
be given. For instance, one may take the Laplace transform of (4.1), (4.2), at 
least formally: the new unknowns satisfy a problem similar to (3.1)-(3.4). This 
means that the interface conditions for the new unknowns are precisely (3.8) or 
(3.9): by anti-transforming these conditions one gets exactly (4.6) or (4.7). 

Furthermore, problem (3.1)-(3.4) can be viewed as a (possible) steady 
state for the time-dependent problem (4.1)-(4.5), or else as the time- 
discretization (at any time level) of problem (4.1 )-(4.5 ), using an implicit 
time-stepping scheme. In both cases, in section 3.1 we have seen that the inter- 
face conditions (3.8) or (3.9) are appropriate for problem (3.1M3.4). Thus, 
(4.6) or (4.7) turn out to be appropriate for problem (4.1)-(4.5). 


4.2. The numerical approximation 

First, we consider a semidiscrete (continuous in time) approximation of 
problem (4.1)-(4.5), endowed either with (4.6) or with (4.7). Keeping the same 
notations of the preceding sections 2.2, 3.2, we apply the spectral collocation 


method in space, that is, we look f or two mappings 

r w^(r ) € (P*) 3 

satisfying, for all t >0 and all j — 1, • • • ,N— 1, 

w ,v.f + A Wn,x + ^ w jv = F at x~ , (4.8) 

— )]* + A w^ >JC + B = F at x + , (4.9) 

At the left boundary we impose the conditions 

G _w jv = q _ , T n [ w^ t + A w^ x + B — F ]= 0 (4.10) 

for x = x q and t > 0, while at the right boundary the conditions are 

G + w jv + H + Wjv * = q + (4.11) 


for x = Xpj and t > 0. Eventually, the two alternative sets of interface condi- 
tions to be requested for x = x^ = Xq and t > 0 are the following: 

(a) Variational approach , 


T p [w N + + A Wjv >x + B — F] = 0, 
Tn w* = T n w^, 


(4.12) 

(4.13) 


— m w jv,x + A w N = A w^y. 


(4.14) 



(b) Nonvariational approach. 


T p [w v , + Aw Nx + Bw n — F] = 0, (4.15) 

w yV — w iV- (4.16) 

Tn Vn,x = T n Wu tX . (4.17) 


A fully discrete approximation to problem (4.1)-(4.5), endowed either 
with (4.6) or with (4.7) can be achieved by applying a time-stepping procedure 
to (4.8), (4.9). Whatever scheme (either implicit or explicit) one uses to advance 
from a known time level t k to a new one t k+1 , the interface conditions, as well 
as the boundary conditions , need to be imposed at the new time t k +1 . 

If an explicit scheme is used in this regard, at the time t k+1 the unknown 
vectors {w^(x^ - )} and {w^U^)}, ; = 1, • • • ,N — 1, can be computed indepen- 
dently of the boundary and interface values. Once these internal values are 
available, the boundary equations (4.10) and (4.11), together with the inter- 
face conditions (4.12)-(4.14) (or (4.15M4.17)), can be solved to provide the 
remaining values at boundary and interface points. Actually, we note that the 
presence of derivatives in space among boundary and interface conditions relates 
boundary and interface values to each other. We also note that the differential 
equations between brackets in (4.10) and in (4.12) (or (4.15)) ought to be 
advanced by the same explicit scheme which was used for the equations at the 
internal points. 

When an implicit time marching scheme is used, the internal unknowns are 
not decoupled f rom the remaining ones any more. As an example, we detail the 
case of the simplest implicit scheme, namely the first order forward Euler 
scheme. 

Denoting by A t the time step, by t k = k At the k -th time level and by 
(w ft ) k the spectral solutions at the time t k , the corresponding problem reads: 

(w A -)* +1 + At [A w^ x + Bwx — F ]* + 1 — (w^)* = 0, (4.18) 

(wjv)* + 1 + At { - *)]* + 

+ ^ + 5 Wjv - F f +1 - (wjf/ = 0. (4.19) 

The boundary equations (4.10) and (4.1 1) are discretized as follows: 


[G w N - q ]* + 1 = 0, 
T n { (w^)* +1 + At [A w„' X + B - F] k+1 - (wjf) k }= 0 


(4.20) 
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[ G + w yy + H + w^ x — q + ]* +1 = 0 (4.21) 

at X\. Analogously, the interface conditions (4.12)-(4.14) give: 

T p \ (w ,y)* + 1 + At [A x + B - FF + 1 - (w^)* }= 0, (4.22) 

T n [(w*)* + 1 - (w^) fc + 1 ] = 0, (4.23) 

— fj. k +1 (w,if,i) H1 + A (w J y)* +1 = A (w^)* + 1 , (4.24) 

at x 0 + . The alternative interface equations (4.15M4.17) read: 

T p \ (w^) ul + At [A w x tX + B — F] fe + 1 — (w^)* }= 0. (4.25) 

(w^)* +1 = (w^)* +1 , (4.26) 

T n [(w.v,J* +1 - (w^ x )* + 1 ] = 0. (4.27) 


We notice that the structure of the system would be the same when using 
other implicit time-marching schemes (such as, for instance, the second order 
Beam & Warming scheme). 

Remark 4.1 We note that (4.18M4.21) with the interface conditions (4.22)- 
(4.24) (or (4.25)-(4.27)) have the same shape as the time independent prob- 
lem (3.33)-(3.42) considered in the previous section. Clearly, in (3.33)- 
(3.42) we must replace w by w jv +1 , B by B + (A?) -1 / and F by 
F* +1 -I- (A? ) -1 wjy, respectively. Therefore, the same iterative procedure can 
be used in order to decouple the hyperbolic problem in Cl~ and the elliptic 
one in fi + . 


APPENDIX: abstract analysis of the regularizing problems presented in 
section 2 


In this Appendix, we detail the existence and asymptotic convergence 
results stated in Propositions 2. 1-2.4 for problems (P € ), (Q € ), (P e )jV anc * (Qe\v • 

As a standard notation, whenever Q is an open interval and k is a positive 
integer we introduce the Sobolev space (see [A]) 

H*(0)= { v € L 2 (Q): Z) m v € L 2 ( Q), m= 1, • • • ,k }. (A.l) 


H* (Q) is a Hilbert space with norm 


H * ( fi) 


k 

L 

m - 1 


lvl& 0 ,+ L II A" V II l H0) 


IV, 


Since Q is one dimensional, we have that 
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H* (ft) C C* _1 (ft ), (A.2) 

for all positive integer k , the embedding being compact. In particular, H K ft) is 
made by functions continuous up to the boundary. Therefore, the following 
(usual) notation is meaningful: 

H 0 1 (fi)= { v € H ! (Q): v = 0 at the endpoints of ft}. (A. 3) 

Unless otherwise stated, we will make the following assumptions on the 
data of problem (P): 

At€L°°(6,c), a€H 1 (a,c), j3 € L 2 (a ,c ), / € L 2 (a ,c ). (A.4) 


Problem (P 6 ). 


Recall that a > 0 in this case. 

Under the assumption (A.4) (actually, under milder assumptions), (P e ) 
can be written in a rigorous variational form: 


find w e € W such that , for all (t>€ W, 


C 


fa e w 


e,x 


c c c 

(t> x dx + Jaw € X (t>dx + J(iw € (l)dx = Jf (t>dx , 

a a a 


(A.5) 


where 


W = Hq 1 (a ,c ), a e = 

If u £ solves (A.5), then the functions 


€ in (a ,b ) 
ft in ( b ,c ). 


U 6 — v € I (a ,6 ) > v € W € I (b ,c ) 


(A.6) 


(A. 7) 


solve (2.5M2.9): this is easily checked by means of suitable choices of 0 in 
(A.5). In particular, (A.5) entails the equation (in the distribution sense) 


— (a € w € x )* + o^e .x + = / in ( a ,c ), (A.8) 

whence 

a e w e x € H 1 (a ,c ). (A.9) 

By (A.2) it follows that both w € and a € w ex are continuous in [a,c], hence 
(2.7), (2.8) and (2.9) have the classical meaning. 

In order to achieve an existence result for (P e ), from now on we make .the 
following requests: 
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M ^ Mo W (* .C ). 

for a suitable strictly positive constant Mo and 

2/3 — a x ^ 0 in (a ,c ). 


(A. 10) 
(A.l 1 ) 


Lemma A.l Under the assumptions (A.4), (A. 10), (A.l 1), (P e ) has a unique solu- 
tion. 

Proof. It is possible to apply Lax-Milgram lemma, because (A. 10), (A.ll) and 
Poincare inequality imply that problem (P € ) is coercive. It goes without saying 
that coerciveness fails as e — ► 0. 

□ 

Now, let us discuss the asymptotic behavior of w € as e — ♦ 0. We recall the 
notations (A. 7) and the assumptions (A.4), (A. 10) and (A.ll), which still hold. 


Lemma A.2 There is a constant C > 0 such that 


IIW € H L 2( a , c ) ^ C ’ 

(A. 12) 

" v e,x Hl 2 ( 6 , c ) ^ ^ ’ 

(A. 13) 

II 1 *L 2 (a ,b ) ^ C- 

(A. 14) 

Proof. Plug the function d> in (A.5), with </> = e~ x 

w € in (a ,b ), d> = e~ b u’ € in 


( h ,c ), then integrate by parts. The assumptions and Poincare inequality give the 
results. 

□ 


Lemma A.3 The L 2 norm of u 6 x is bounded in a right neighborhood of the left 
boundary x = a . 

Proof. Let xlt be a smooth function in (a ,c ), vanishing outside a right neighbor- 
hood of a. Take the L 2 (a ,c ) scalar product of (A.8) by \lrw ex : the assertion 
follows by (A. 12) and (A. 14). 

□ 


Now, let us introduce the function 

= a t w e,x 

We already know that $ e € H Ka ,c ) (see ( 


— aw ( : 


(A. 15) 
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Lemma A.4 The H 1 norm of <i> € is bounded in (a ,c ). 

Proof. Lemma A.2 gives the boundedness of <& € in L 2 (a ,c ); (A.8) and (A. 12) 
give the boundedness of <I> € x in L 2 (a ,c ). 

□ 

Now, we are in a position to give the following result, which completes 
and refines the statement of Proposition 2.1. 


Proposition A.1 Assume (A.4), (A. 10), 
v € L 2 (b ,c ) which satisfy 

(A.ll). There are 

u € L 2 (a ,b ) and 

a u x + (3 u = f 

in L 2 (a ,b ); 

(A.16) 

~ (mv x ) x +(*v x +£v = / 

in L Kb ,c ); 

(A.17) 

K 

cT 

II 

o 


(A. 18) 

v (c ) = 0; 


(A.19) 

OiU = — yUV x + QfV 

at x — b. 

(A.20) 


Proof. As a consequence of Lemmas A.2-A.4, of Banach-Alaoglu-Bourbaki 
theorem and of (A.2), we can find u€L 2 (a,h), v€L 2 (6,c) and 4>€H 1 (a,c) 
such that (upon extracting a subfamily) 


(i) 

u e — 

u 

weakly in L 2 (a ,b ); 

(ii) 

V € —* 

V 

weakly in H Hb ,c ); 

(iii) 


<t> 

weakly in H x (a ,c ); 

(iv) 

u e (a ) 

— » 

u (a ); 

(v) 

v € (c ) 

- 

v(c ); 

(vi) 

v e (b) 

- 

v(fi ); 

(vii) 

&e,x 

-» 

0 strongly in L 2 (a ,b ). 


Note that the value u (a ) is well defined, because of Lemma A. 3. (i)-(iii) and 
(vii) permit to pass to the limit in (A.5): this gives (A.16) and (A.17). (A. 18) 
and (A.19) follow by (iv) and (v), respectively, since u € (a ) = v e (c )= 0. 
Finally, (i)-(iii) and (vii) entail that O = — ocu in (a ,b ) and $ = fiv x — av in 
(b ,c ), whence (A. 20) follows, by (vi). 


Remark A.1 Analogous results could be proved when replacing the 
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homogeneous Dirichlet condition at c by a Neumann condition or by a 
Newton-type condition. 


Problem (P € ) A r 


Recall that a < 0 in this case. 

For this problem, the variational formulation is still (A.5), just changing 
the function space: now we take 

W = { v £ H ’(a ,c ): v (c ) = 0 }. (A.21) 

The existence holds the same way as in the previous case and the asymptotic 
analysis is analogous. We detail the main steps, under the assumption: 

0€L “(a ,6). (A.22) 

Moreover, we still assume (A.4), (A.10), (A. 11) and use the notations (A. 7). 


Lemma A.5 There is a constant C > 0 such that 


11 W € l! L 2 (a ,c ) ^ C ’ 

(A.23) 

v «,x ^L 2 (d> ,c ) ^ C ’ 

(A.24) 

II u t,x Hl 2 (<2 ,b ) ^ C ’ 

(A.25) 

V7lu € , x (6)l <C, 

(A.26) 

(Mv e>x ) x ll L 2 (6 c) < c. 

(A.27) 


Proof. (A.23) and (A.24) follow by plugging the function <t> in (A.5), with 
(f)=e x M e in {a ,b ), <f)=e b -w e in (b ,c ), then integrating by parts. The 
assumptions and Poincare inequality give the results. 

To prove (A.25) and (A.26), take the L 2 (a ,c ) scalar product of (A.8) by * , 
where xlr is any smooth function, vanishing in (h ,c ), then integrate by parts. 
Finally, (A.27) follows by (A.8), (A.23) and (A.24). 


□ 


Thus, we are in a position to prove the main result, which was summarized in 
Proposition 2.2. 


Proposition A.2 Assume (A.4), (A. 10), (A.ll), (A.22). Moreover, assume that /j. 
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is continuous at x = b . Then, there are u € H Ha ,b ) and v £ H \b ,c ) which 
satisfy (A. 16), (A. 17), (A. 19) and the interface conditions 

u(b) = v(b), (A. 28) 

v x (b) = 0. (A. 29) 


Proof. As a consequence of the previous Lemma and of (A.2), we can show the 
existence of u € H Ha ,b ) and v € H Hb ,c ) such that (upon extracting a subfam- 
ily) 


(i) 

u t — 

u weakly in H Ha ,b ); 

(ii) 

v e — 

v weakly in H 1 {b ,c ); 

(iii) 

P V €,X 

— > /jlv x weakly in H Ha ,c i 

(iv) 

V £ (c ) 

- v (c ); 

(v) 

u e (b ) 

— > u (b ) and v e (6 ) -+ v(b 

(vi) 


)(b ) -» (/iv x )(b ); 

(vii) 

tUe.x & ) — ► 0. 


(i)-(iii) permit to pass to the limit in (A.5). The conditions at x = c and x ~ b 
follow by (iv)-(vii), noting that ) > 0 (see (A. 10)). 


Remark A.2 If we take a homogeneous Dirichlet condition at x = a instead of 
the Neumann one, then (P e ) A : coincides with (P £ ); so does its variational 
formulation. But now we are assuming a < 0, hence the asymptotic 
behavior is different from that of the case a > 0. It is easy to see that the 
final Proposition A.2 still holds, with u found in L 2 (a ,b ): actually, the con- 
vergence of u € to u is only L 2 (a ,b ) (weak), whence we cannot have a con- 
vergence of u £ (a ) to u (a ), in general. Actually, Figure 2.5 shows a numeri- 
cal evidence of a boundary layer for u e at x = a , although the limit func- 
tion u is obviously continuous in [a ,b] (see (A. 16) and (A.2)). This feature 
makes (P £ ) A - preferable, especially in view of the applications to systems 
(sections 3 and 4). 


Problems ( Q £ ) and (Qe)w • 

Now, the two problems do not admit a "natural” global variational formu- 
lation and the question of existence and the asymptotic behavior are somewhat 
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more complicate. Nevertheless if we assume that (A.4) holds and that 

/j. is continuous at x = b , (A. 30) 

then the equations and the boundary and interface conditions defining (Q e ) and 

(Q € ) a - make sense, provided the solutions are sought for in HKa.b) and 
H Kb ,c ), respectively. 

We begin with problem (Q £ ), recalling that a > 0. It can be shown that 
(O e ) has a unique solution u e , v e , at least for e small, under the assumptions 
(A.4), (A. 10), (A.ll), (A.30). The asymptotic behavior is being investigated 
now, under the same assumptions. 

Lemma A.6 There is a constant C > 0 such that 

HMlW) ^ C, (A. 31) 

(A - 32) 

lw e (b)l<C. (A. 33) 

Proof, (i) Take the L 2 (a ,b ) scalar product of (2.5) by n{b)e~ x u e , then 
integrate by parts. 

(ii) Take the L 2 (b,c) scalar product of (2.6) by ee~ b v e , then integrate by 
parts. 

(iii) Add the two equations provided by (i) and (ii), term by term: the conclu- 
sion follows by Poincare inequality. 

□ 


Lemma A.7 There is a constant C > 0 such that 


** V e ^H'(6,c) ^ C • 

(A. 34) 

1 v e x (b ) 1 < C, 

' (A. 35) 

" '(a ,b ) ^ 

(A. 36) 


Proof. Let £ £ € H Kb ,c ) be the solution of 

_ (M€.;c)x = 0 in (b,c), £ € (b )= v £ (b ), £ £ (c)=0. 

By (A.33), the HKb ,c ) norm of £ £ is bounded, as well as the value of 4 e>;c (b ). 
Moreover, the function d £ = v £ — belongs to H 0 Kb ,c ) and satisfies 


)x + ocd £Jt + fid € — g e , 


(A. 37) 
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where g e = f — ul, ex — /K« is bounded in L 2 (h ,c ). Multiplying (A. 37) in 
L 2 (b ,c ) by d e , it follows that the H 1 (b,c) norm of d e is bounded, whence 
(A. 34). 

Next, we multiply (A. 37) by \lf/xd e x , where xls is a smooth function vanishing 
outside a right neighborhood of b : (A. 35) follows easily. 

Finally, (A. 36) can be proved by taking the L 2 (a ,b ) scalar product of (2.5) by 
x and using (A. 31), (A. 35). 

□ 


From Lemma A. 7 we get the following proposition (see Proposition 2.3). 

Proposition A.3 Assume ( A.4 ), (A. 10), (A.ll), (A.30). There are u £ H*(a,h) 
and v € H x (b ,c ) which satisfy (A. 16), (A. 17), (A.18), (A. 19) and (A.28). 

Proof. Let u e , v e solve (Q e ). By Lemma A. 7, there are u€H 1 (a,h) and 
v € H 1 (b ,c ) such that (upon extracting a subfamily) 

(i) u e — ► u weakly in H Ha ,b ); 

(ii) v e — ► v weakly in H Hh ,c ); 

(iii) u e (a ) -♦ u (a ); 

(iv) v 6 (c ) -► v(c ); 

(v) u e (b ) — > u(b ) and v e (b ) -♦ v(b ). 

All of these properties permit to pass to the limit in the regularized problem 
(Q e ). Thus, the proof follows easily. 

□ 


Now we come to problem (Oe^ : recall that a < 0. 

This case looks somewhat trickier than the previous one and the natural 
choices for test functions do not seem to be appropriate, in proving the a priori 
estimates. Even more, it can be shown that problem (Oc)a' ma Y fail t0 have a 
solution under the assumptions (A.4), (A.10), (A.ll), (A.30) (which were 
sufficient f or existence in the previous case). 

This trouble seems to be motivated by the lack of a maximum principle 
under the sole coerciveness condition (A.ll) on 0. For this reason, we discuss 
problem (Qe)^ under the further hypothesis: 

0(x ) > 0 for x a.e. in (a ,b ). (A. 38) 

We just note that such an assumption is not strongly restrictive if the problem 
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we are dealing with is regarded as a time discretization of an evolution problem 
by an implicit method (see section 4.2). 

It can be shown that (Qe).v has a unique solution u e . v e , under the 
assumptions (A. 4), (A. 10), (A. 11), (A. 30), (A. 38). 

The asymptotic behavior is being investigated now, under the same 
assumptions: for technical reasons, we will confine the situation a bit more, 
making the further hypothesis: 

/ £ L°°(a ,6 ), /3(jt)^/3 o >0 for x a.e. in (a ,b ), (A. 39) 

for some /3 0 . This allows us to get low order estimates on u e and v € . Later on, 
we will make f urther assumptions in order to find higher order estimates. 

Lemma A.8 There is a constant C > 0 such that 


11 V € ll H 1 (6,c) ^ C ’ 

(A.40) 

\v ex (b)\ < C, 

(A.41) 

) x N L 2 (6 >c) < C. 

(A.42) 


Proof, (i) Take the L 2 (a ,b ) scalar product of (2.5) by nib )<£u e , where 

6 

d)(x ) = — exp {— fa(t )dt } . 

Then, integrate by parts. 

(ii) Take the L 2 {b ,c ) scalar product of (2.6) by v e , then integrate by parts. 

(iii) Add the two equations provided by (i) and (ii), term by term. Recalling 
(A. 10), (A.l 1 ), (A. 39), we find that 

be b 

e/xib ) fd>u e x 2 dx + yU 0 Jv e x 2 dx + fiib ) $ 0 J(bu 2 dx — y aib )v 2 {b ) < 

a b a * 

b c 

< /xib ) Jf (t>u e dx + £fv € dx . (A.43) 

Now, by (A. 39) we have 

b b 

Jf <t>u e dx < 11/ ll L ~( ? l“e }dx , 

a a 

so that Poincare inequality in (A.43) gives 

b c b 

e fd>u ex z dx + k ifv e X 2 dx + ^ 0 f(t>u e 2 dx +k 2 v e 2 ib )< 

a b a 
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\u e \dx + k 3 , 


(A. 44) 


where k t are positive constants, i = 1,2,3. In particular, it follows that 

b b 

( 3 0 J<t>u e 2 dx ^ II / H L ~( a lu £ IJ.X + £3 

a ’ a 


b 

and an elementary computation shows that the integral Jd> \u € \dx is bounded. 

a 

Thus, (A.44) and Poincare inequality imply (A. 40) and the boundedness of 
v £ (6 ). 

To show (A.41), take the L 2 (b,c) scalar product of (2.6) by fi\lfv € x , 
where x[r is a nonnegative, smooth function, vanishing near c, with t|f( 6 )= 1 . 
After integration by parts, (A.41) follows by (A.40). 

Finally, (A.42) follows by (A.40) and by the very equation (2.6). 

□ 


Lemma A.9 There is a constant C > 0 such that 

ll ^ ll H'(a, t ,)< C * (A ’ 45) 

Proof. Take the L 2 {a,b) scalar product of (2.5) by e x u e , then integrate by 
parts. By Lemma A. 8 -, it follows that the L 2 (a ,b ) norm of u e is bounded, as 
well as the value u € (a ). 

Next, take the L 2 (a ,b ) scalar product of (2.5) by u € x , then integrate by 
parts. The conclusion follows by Lemma A. 8 and by the first part of this 
proof. 

□ 

Now, we are in a position to prove part of the results stated in Proposition 

2.4. 

Proposition A.4 Assume (A.4), (A. 10), (A.ll), ( A.30 ), (A.39). There are 
u € H *(a ,b ) and v € H l (b ,c ) which satisfy (A. 16), (A.17), (A. 19), (A.28). 

Proof. Let u e , v € solve (Q e ) N • By Lemmas A.8, A.9, there are u € H Ha ,b ) and 
v € H Kb ,c ) such that (upon extracting a subfamily) 

(i) u e —>u weakly in H ! (a ,b ); 
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(ii) v € — * v weakly in H ,c ); 

(iii) v € (c ) — > v (c ); 

(iv) u e (b ) — > u (b ) and v e (b ) -» v (b ). 

All of these properties permit to pass to the limit in the regularized problem 
(Q e ) A ;. Thus, the proof follows easily. 

□ 


To complete Proposition 2.4, it remains to show that the derivatives of the limit 
functions u,v of the preceding Proposition join continuously. To this end, we 
assume that the data a, (3, f are more regular than it was until now, precisely: 


a is Lipschitz continuous in [a ,b ], 


0€H \a,b), f €H \a,b). 


(A. 46) 


Lemma A.10 Assume ( A.4 ); (A. 10), (A.ll), ( A.30 ), (A.39), (A.46). The L 2 norm 
of u e xx is bounded in a left neighborhood of the interface point x = b . 

Proof. Take the L 2 (a ,b ) scalar product of (2.5) by u e xx (which lies in L 2 (a ,b ) 
because of the equation itself). Next, integrate by parts in all terms except in 
the first. Recalling (A.41), (A.45) and (A.46), we get 

(A>47) 

for some C > 0. Finally, take the derivative of (2.5) and multiply it in L 2 (a ,b ) 
by <t>u €xx , where 0 is smooth, nonnegative, with 0(a )= 0. By (A.45), (A.47) 
and recalling that a < 0, the assertion follows. 

□ 


Proposition A.5 Assume (A.4), (A. 10), (A.ll), (A.30), (A.39), (A.46). The func- 
tions u , v considered in Proposition A.4 satisfy 

u x = v x at x = b. (A.48) 

Proof. Since the property holds for u e and v € (see (2.12)), it is enough to prove 
that: 

(i) u ex (b ) -> u x (b); 

(ii) v € x (b ) — » v x (b ). 

(i) follows by Lemma A. 10 and by (A.2); (ii) follows by (A. 2) and by (A.42), 
recalling (A. 10) and (A.30). 

□ 
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Thus, the proof of Propositions 2. 1-2.4 is complete. 
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